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NATIOML ADVISORY CCMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1291 


ON STABILITY MD TURBULENCE OF FLUID FLOWS* 
By Werner Heisenberg 


INTRODUCTION 


The t\irbulence problem, which will quite generally form the subject 
of the following investigations, has been treated in the course of time 
in so many reports from so many different viewpoints that it is not oiir 
intention to give, as an introduction, a survey of the results obtained 
so far. For that piarpose, we refer the reader to a report by Noether^ 
on the present state of the turbulence problem, where most biblio- 
graphical data may be found as well. 

For our purpose, a rough outline of the present state of the tur- 
bulence problem will be svLfficient. The investigations made so far are 
divided into two parts; one part deals with the stability investiga- 
tion of any laminar motion, the other with the turbulent motion itself. 

The first-mentioned investigations led, at the beginning, to the 
negative res\ilt that all laminar motions investigated are stable. 

V. Mises2 and L. Hopf3 proved, on the basis of a formula by Sommerfeld^, 
the stability of the linear velocity profile corresponding to Couette*s 
arrangement. Blumenthal^ reached the same result for a profile of the 
third degree, upon which Noether invited discussion. On the other hand, 
Noether^ later succeeded in specifying an unstable profile - a profile 
which is unstable even in the case of a frictionless fluid can never be 
realized as a steady state of motion for actual conditions. More 


*'*Uber Stabilit'dt und Turbulenz von Flussigkeitsstrdmen. " Annalen 
der Physik, Band 7^, No. 15, 192^, pp. 577-627. 

^Noether, F. : Zeitschr. f, angew. Math. u. Mech. 1, 1921, p. 125. 

%ises, R. V.: Beitrag z. Oszillationsprobl. : Heinr. Weber- 

Festchrift, 1912, p. 252. 

^Hopf, L. : Ann. d. Phys. 44, 1914, p. 1. 

^Sommerfeld, A.; Atti d. IV. congr. int. del Mathem. Rom 1909. 
^Blumenthal, 0.; Sitzxmgsber. d. bayr. Akad. d. Wiss., 1913^ P. 563. 
Noether, F. : Nachr. d. Ges. d. Wiss., Gottingen 1917. 
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recently, however, Prandtl"^ has shown that indeed profiles exist which 
possess unstable characteristics only if the friction is taken into 
consideration. 

The other group of reports which achieved great success quite 
recently by the calculations® of Von Kerman, Latzko, and others 
investigates the turbulent motion itself proceeding by a semiempirical 
method using the laws of similarity. Theoretically, the reports of 
this group are based almost throughout on Prandtl’s boundary-layer 
theory. Their most important result for our purpose is the so-called 
yl/7 -law of turbulent velocity distribution which follows from Blasius' s 
law of resistance (examples can be found in Schiller's report^.) 

The determination of the critical Reynolds number was always one 
of the main aims of the first-mentioned reports, the stability inves- 
tigations. So far, a satisfactory calculation of this number has not 
been accomplished and it must be regarded as doubtful whether it could 
be achieved by stability investigations. The tests of Ekman^O, RuckesH^ 
and Schiller^, together with the negative results of Hopf concerning 
the linear velocity profile, rather suggest the notion that the critical 
Reynolds number does not indicate the point where the laminar motion 
starts to become unstable, but the point where, for the first time, the 
turbulent motion is possible as steady state. From the viewpoint of 
theory, we must thus be prepared to find eventually two critical 
Reynolds numbers, one corresponding to the beginning of turbulence, the 
other to the breaking down of the laminar motion. 

The present investigation also will be divided into two different 
parts, the treatment of the stability problem on the one hand, that of 
the turbulent motion on the other. 


Tprandtl, L. : Physik. Zeitschr. 23, 1922, p. 19, and Tietjens, 0., 

Dissert. Gottingen, 1922. 

®See Zeitschr. f. angew. Math. u. Mech. , 1, 1921, p. 233f* 

9schiller, L. : Rauhigkeit \md"^ kritische Zahl. Physik. Zeitschr. 3^ 

1920, p. kl2. 

^O pknign , V.: Turbulent Motions of Liquids. Arch. f. Mat. och 

fysik 6, 1919* P* 12. 

llRuckes, W. : Dissert. ‘Wurzburg I9OT. See also Lecture by W. Wien, 

liber turbul. Bewegungen. Phys. Zeitschr. 8., 190^; and Verh. d. deutsch. 
phys. Gesellsch. 9> 1907 • 
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The aim of the first part is to summarize all previous investiga- 
tions under a unified point of view, that is, to set up as generally as 
possible the conditions under which a profile possesses tinstable or 
stable characteristics, and to indicate the methods for solution of the 
stability equation for any arbitrary velocity profile and for calcula- 
tion of the critical Reynolds number for unstable profiles. This aim 
can, of course, be attained only imperfectly by the use of approxima- 
tion methods. Nevertheless, we hope to be able to clarify by such 
calculations the qualitatively essential viewpoints. At first, the 
investigation of any arbitrary profile seems physically meaningless 
since only certain profiles actually occur; however, since we may 
interpret any profile as finite disturbance of another, as for instance 
Noether has done elsewhere, and since we must, on the other hand, later 
extend the investigations to the (at first unknown) basic profile of the 
turbulent motion, the investigation of an arbitrary profile seems, after 
all, to be of great importance. 

As application of the methods, the parabola profile will be cal- 
culated completely. 

In the second part, we shall attempt to derive, tinder certain 
greatly idealizing assumptions, differential equations for the turbulent 
motions and to obtain from them qualitative information about several 
properties of the turbulent velocity distribution. 


PART I: THE STABILITY EQUATION 

1 . Statement of the Mathematical Problem 


The most essential limitation we impose on our calculations con- 
sists in the exclusive consideration of two-dimensional laminar motions 
and only two-dimensional disturbances of these motions. Taking a rec- 
tangular coordinate system X, Y, Z as basis, we therefore assume that 
the velocity in the Z direction is zero and all remaining quantities 
independent of Z. Furthermore, however, we shall only examine the sta- 
bility of such laminar motions as occvir between two straight parallel 
walls. We assume the walls to be parallel to the X axis; therefore, 
the laminar motion to be investigated also promises a velocity com- 
ponent only in the X direction. This velocity w in the X direction 
will, in some way, be dependent on y. Concerning the function 
w = w(y) we reserve for later making a few assumptions about continuity, 
symmetry, etc.; otherwise, however, this function is to be at first 
quite arbitrary. 
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If we put w = ay, our formulations become exactly identical with 
those investigated by Hopf in the Couette case. The problem whether 
the investigated profiles w = w(y) can be realized as steady motions 
will not be dealt with for the present. 

Before deriving once more the stability equation (already set up 
elsewhere by Sommerfeld) briefly from Stokes’s differential equations, 
we introduce dimensionless variables in the known manner. Let h be 
a characteristic length (for instance the distance between the two 
walls), U a characteristic velocity of the profile, p, the viscosity, 

p the density, and — ^ = R the Reynolds number; we introduce instead 

P 

of X, y, u, V, t, and p (u, v being the velocity in x or 
y direction, respectively, t the time, and p the pressiire) new 
variables x^, y^, u^, v^, t^, and p^, according to the relations 
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By A one understands here the differentiation symbol 


Equation (^i-) does not yet contain anything about our special 
problem, the stability investigation of a certain laminar flow. Accord- 
ingly, equation (4) will form also the basis for the calculations of 
part II. In order to pass over specifically to the stability investi- 
gation, we divide the motion and therewith also the vector potential i|r 
into a basic flow and small oscillations superimposed over it. Thus we 
set up the formula 


3>(y) + cp(y)ei(P^-<^) 

(5) 

^ = w(y) = w 

oy 

(6) 


If we enter this formula into equation (4), omitting all terms not con- 
taining cp (since we regard equation (4) as satisfied for <p = O), 
furthermore omitting all terms quadratic in cp ( since we assume cp as 
small), the corresponding differential equation for q> reads 


(cp'' - a^cp)^w - “ Tw" " ^ ■ 2a^cp" + a^) (T) 


The fact that we regard equation (4) as satisfied for cp = 0 signifies 
physically that we consider only such basic flows w which either, by 
virtue of external forces, are really steady, or show a variation with 
time which is slow compared to that of the small oscillations. 

Equation (?) is in this generality already derived elsewhere by 
Noether. It is an ordinary differential equation for cp of the fourth 
order. It corresponds to the fact that the function cp must fulfill 
four boundary conditions; u and v , thus also (p and cp* , must dis- 
appear at the two walls. If we put p/a = c so that c essentially 
signifies the wave velocity, the mathematical problem may be formulated 
as follows: The solutions of the equation 
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(cp" - a^cp)(w - c) - cpw” = — (cp"" - 2a^cp" + a^cp) (Ta) 

aR 

are to be investigated with the secondary condition that at the bovmding 
walls (for instance, y = 1 and y = -l) <p = 0 and 9 ' = 0. For each 
value of a and R the corresponding value of c and 3 is to be 
calculated; let a for reasons of simplicity always be positive. 
According to whether the imaginary part of 3 is positive, zero, or 
negative, we are dealing with a stable, undamped (undamped oscilla- 
tions = neutrally stable oscillations), or unstable oscillation. The 
conditions for profile w are to be fo\ind under which equation (Ta) 
admits only stable oscillations or, respectively, also unstable ones. 

Before turning to the methods of solution we want to point out a 
special property of the equation (Ta). In the limit of the frictionless 
fluid, R = 00 , equation (Ta) is transformed into a differential equa- 
tion of the second order for cp 


(cp" - a^cp) ( w - c ) - cpw" = 0 


(8) 


Accordingly, only two boundary conditions must now be satisfied which 
signify that the normal velocity component, thus v or cp but no 
longer cp' , is to disappear at the two walls. 

The conditions for the solvability of equation (8) have already 
been investigated in detail by Rayleigh. Introducing a simple desig- 
nation, one may distinguish basic flows "capable of oscillation" or 
"not capable of oscillation" according to whether or not equation (8) 
possesses a solution with real c which satisfies the boundary con- 
ditions. jf solutions with complex c exist, the stability problem 
for these oscillations has, as will be shown later, already been 
decided by equation (8), also in case of consideration of the friction; 
the oscillations are then always unstable. 

One is, however, beyond this led to the conjecture that the pro- 
file w, under influence of friction, permits imstable or inadamped 

^^Lord Rayleigh, Papers I., p. 36 I; III. pp. 5T5^ 59^; IV. p. 203. 

^^Here, however, it is by no means sufficient to approximabe the 
profile by tangents polygons; the result with respect to possible 
oscillations would thereby be completely falsified. 
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oscillations only in one case: when it "belongs to the "basic flows 

capable of oscillation. 

This supposition is the more obvious as it has been confirmed for 
all profiles investigated so far.l^ Nevertheless, it is by no means 
motivated by the fact that equation (8) results from equation (7a) in 
the limit R = <» since it has been proved, for instance, in the reports 
of 0seen^5 that the limiting process R = » has led more than once to 
false results in the differential equations, particularly with respect 
to the boundary conditions of the frictionless fluid, and that one may 
therefore apply the limiting process only to the integrals of equa- 
tion (7a). Moreover, it can by no means be decided beforehand whether 
the friction modifies the imdamped oscillations of equation (8) in the 
sense of a damping or an excitation. 

Following, we shall attempt to prove our surmise mentioned above by 
showing that the systems capable of oscillation are shown to possess 
above a certain value of the Reynolds number and in general unstable 
character, whereas all systems not capable of oscillation are shown to 
possess a stable character. 

By this principle the problem of the stability of a profile is 
quite considerably simplified since, as is well known, the solutions of 
equation (8) may be directly written down for very small values of a. 


2. The Methods of Solution and the General Behavior of 
the Integrals of Equation (7a) 

The most important property of equation (7a) which permits an 
approximate representation of its solutions consists in the fact that 
R may be regarded as very large. It will become evident that if a 
stability limit exists, this limit lies, in general, at very high values 
of R. Since it is, moreover, physically quite improbable that for 
small values of R instability of the respective profile could occur, 
it is siafficient for our next purpose to regard R as very large. 

This assumption makes it possible to approximate the solutions of 
equation (7a) by development in negative powers of R or, as will be 

shown, of Furthermore, we shall assume a as small and shall 

develop the solutions in a given case in positive powers of a^. 


^^Compare also Prandtl, Physikalische Zeitschrift. 

^^Compare, for instance, C. W. Oseen, Beitrage z. Hydr. Annalen 
der Physik k6, pp. 231 and 623# 1915* 
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Both methods of development in and respectively, 

seem contradictory insofar as in the first case oR is assumed large, 
in the second case small; however, the contradiction is eliminated 

by the fact that R may be regarded, in general, as extraordinarily 
large so that for instance for R = 2000, a = 1/10, aR becomes equal 
to 200, a? = 1/100 which is fully sufficient for a satisfactory con- 
vergence of the two developments. However, the convergence properties 
of these approximation methods must be considered more exactly. The 

investigation shows that the series in are generally diver- 

gent, yet show the well-known characteristics of the semiconvergent 
series, that is, that the terms first decrease, then again increase, 
and that one obtains the optimum approximation if one breaks off the 
series with the smallest term. Our approx ima.t ion method has, there- 
fore, convergence properties similar to those of the series of the 
perturbation theory used in astronomy, the behavior of which is 
described in detail by Poincar^, Meth. nouv. d. 1. mec. cel. II. 

The use of the semiconvergent developments is rendered con- 
siderably difficult by the fact that they lose their validity in the 
neighborhood of a certain point so that it cannot be immediately 
decided in what manner the approximate solutions on both sides of the 
point must be joined in order to approximate a certain integral of the 
equation (Ta) on both sides. This question will be discussed in detail 
in section 3. 

The development in positive power series of seems, in general, 

to be actually convergent. For special profiles this development may 
be strictly proved (for instance, for the linear profile); however, we 
have not carried out an investigation of the problem under what con- 
ditions for the profile w this convergence actually occurs. 

We start with the derivation of the approximate solutions of 
equation (7a) 


(cp" - a^cp)(w - c) - w''cp = ■ Sa^cp" + a^cp) (7a) 


For this purpose we first put 


/gdy 


y 


g = \/aBgQ + g^ 



cp = e 


2 




(9) 
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We shall limit the development to the two highest terms in \JaR. There 
follows 


oKgO^(’W - 



+ 2gogi)(w - c) = ioRgQ^ + lySS 


and w" are presupposed to be of the order of magnitude 1 or, at 
any rate, « ^"oR. By means of simple calculation there now results 


Sq = - c), g^ = - I dy = - I log gQ 


( 10 ) 


Thus we obtain two particular integrals of the equation (7a) 


/, * / )/^i(xR(w-c) 

. (V . cj-SAeJl'o' 


( 11 ) 


The point is to be determined by w being = c for y = y^. 

Thus Yq may, under certain conditions, be complex. The sign of the 
root is to be chosen so that for 


V - c 


-ae 


ixD 


-i(xR(v - c) = oRae 



the root becomes 




A remarkable fact about these two integrals is that in qp 

does not appear in this approximation (that is, only in the combina- 
tion (xR which, as follows from eqxiation (?a), could in a certain 
sense, be called the true Reynolds nnmber). 

As we shall see later, the two integrals (ll) determine the 
behavior of cp in the boundary layer, and the nonoccurrence of 
in equation (ll) signifies physically that we consider only oscillations 
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the wave length of which is large compared to the boundary-layer thick- 
ness^ which is certainly the case for the empirically observed unstable 
oscillations. 

For a complete system of solution of equation (7a)^ however^ we 
need two further integrals; naturally^ we shall choose those which 
result from the integrals of equation (8) by the development of 

(clR)"^ in a power series. 

For this purpose we first solve equation (8) by the development 
of aP- in a power series. Thus we put 

(cp” - a^cp)(w - c) - cpw" = 0 (8) 

cp = + a^cp(^) + + . . . ( 12 ) 

Hence follows 

cp(0)"(w - c) - q)(®)w” = 0 
cp(l)”(w - c) - cp^^^w” = cp^^)(w - c) 

(p(2) (v - c) - q>^^^w" = cp^^^(w - c) 

By the method of variation of the constants the two integrals 
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result. In addition, these integrals have now to be corrected by 

quantities of the order (<xR)“^, . . . etc., if they are to satisfy 
equation ( 7 a). 

Without writing the corresponding series development down in 
detail, we give as result cp with the quantities of the order (aR)-l 



ctRu (w - c)^ dy^ 


c) + 



qP 4 = (w - c) 






(li^) 


With equations (ll) and (l 4 ), a complete solution system of the equa- 
tion (7a) has been found approximately. 

Before applying this system of solution for satisfying the boun- 
dary conditions, it will be useful to clarify the physical significance 
of the four integrals cp^, cp^, cp^^ ‘to anticipate a few 

results which we cannot establish until later. 

The integrals cp^, cpg are very rapidly variable for the high 
values of R which are of interest to us, as can be seen from the 
exponent of the order \/^. If, therefore, for instance cp^^ is at 

one wall of the order of magnitude 1, it will vanish exponentially at 
some distance from the wall, (in itself, it also could become extraor- 
dinarily large; however, this is naturally prevented by the boxindary 
conditions.) Consequently, cp is composed, except for the immediate 
neighborhood of the wall, of CP3 and cpi^ alone, that is, is 
very similar to the behavior of cp in the frictionless fluid. 
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The fact that does not explicitly occur in cp]_, (com- 
pare equation (ll)), hut does appear in 94 (compare equation l4) ) 

must evidently he interpreted physically to signify that, if is 

assumed to he not very large (a^ < \f^, compare equations (7a.) and (9))^ 
the wave length may he considered as infinitely large compared to the 
houndary-layer thickness, hut not compared to the width of the channel. 
The characteristic difference between the "ho\mdary- layer integrals" 
cpp, cp 2 ; hand, and the integrals corresponding to the friction- 

less fluid cp^, 94 , on the other, is therefore significantly expressed 
in the occurrence or nonoccurrence of a^. 

Concerning the convergence of the development in power series 
of a^, we may hope that for values of of the order of magnitude 1 

it is still amply s\ifficient to enable a good approximation, since for 
a linear profile the series for qp^.^ "Ph tiecome series of the type of 

power development of cos a which in the neighborhood a = 1 still 
converges very rapidly. 

The flow pattern to he expected after all these conclusions corre- 
sponds to the formulations made in Prandtl' s boundary- layer theory. 
Except for the Immediate neighborhood of the walls, the motion obeys 
very nearly the differential equation of the frictionless fluids. To 
the walls themselves, however, adheres a boundary layer the thickness 

of which is of the order of magnitude (<xR)"^/^. In this boundary 
layer the velocity u decreases toward the wall rapidly toward zero 
whereas v is almost zero even outside of the bovindary layer. 


3 . The Connecting Substitutions 

If we want to study the course of the integrals of equation (Ta) 
from one bounding wall to the other, we must take into account that, 
at a point y = Jq in the channel w - c = 0 (or that at least the 

the real part of w - c is zero), therefore the wave velocity there 
agrees with the velocity of the basic flow. At this point, the 
approximation formulas (ll) and (l4) for the integrals of equation (7a) 
cease to be valid. 

Thus it is necessary to know the connecting substitutions for Cfq^, 
cpg, 9^^ and 9j^ which have to be applied for the transition from 

Re(w - c) > 0 to Re(w - c) <0. For this purpose, we develop w 
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and cp in the neighborhood of the critical point Jq in the power 
series of (oR )“^/3 and put therefore 


y - yo 


tj(oR) “^/3 


Furthermore we assume the imaginary part of to be of smaller 

order of magnitude than (aR)“^/^. If it is of higher order of magni- 
tude, the connecting substititions are self evident because then 
nowhere in the entire range of real y does a. "critical point" appear. 
If the imaginary part is of the same order of magnitude, the behavior 
of cp and w may be easily interpolated from the two limiting cases 
just mentioned. At first we may even put 


Im(yo) = 0 


since cp in our case 


Im(yo) « 

may be developed in power series of Im^yQ^ and at first only the 
behavior of 9 for Im^yQ^ =0 is needed. 

Thus we now put 


cp = (Pq + ecpi + e q?2 + 


e = (oR)"V3 


2 2 

w - c = ear; + c brj + . . 



NACA TM 1291 


Then there results from equation (7a.) 

Cpo"" + €cpi"'' + . . . = ecpQ"bq^ - 

2c(Pq'^ + . . . 

Thus in first approximation 


^ tut . ^ _ 

9q ” 


(15) 


in second approximation 

cp^"" = + q>Q"bq^ - 2^^^ (l5a) 

For the integrals cp^, equation (ll), we infer from equa- 

tion (15) that they heha.ve in the critical range (q order of ma,gni- 
tude 1) like the integrals found by Hopf for the linear profile, that 
is, like certain cylindrical functions. Thus, we may conclude at 
this point that the connecting substitutions for cp^, cp^ from equa- 
tion (11) except for quantities of the order (aR)“l/3 must be the 
same as for the linear profile 


Tl — 

9l - i<P2 — > 


(16) 


corresponding to a transition of 


Re(w - c) < 0 ^ Re(w - c) > 0 
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However, for the study of the integrals cp^, in the neighbor- 

hood of y - Yq = 0, the simple calculations made so far are not suf- 
ficient since for the latter the approximate solution (15) would 
read qp" = 0; however, we know from equation (l^<-) that, in the limit 
R — cpij,’ , in general, becomes logarithmically infinite at the 

point y - Yq = 0. Equations (15) and (l5a) are therefore in this 
form unsuitable for expressing this singularity. 

Instead, we now set a = 0 and w"* =0 (that is, we break off 
the development of w with the second term); otherwise, however, 
integrate equation (7a) exactly at first. In doing so, we notice that 
cp = w - c must be a particular integral of this simplified equa- 
tion (7a) and we make, therefore, for qp the statement familiar from 
the theory of linear differential equations 


9 


(w - c) 



i dy 


Then there follows from 


cp"" = -ioR(cp"(w - c) - w"qp) 


for 


cp=(w-c)J^ i|fdy 

■f"' (w - c) + 4V'w' + 6\lf'w" = -ioR(2w' (w - c)t1c + (w - c)^i|r') 


which after repeated integration becomes 


V(w - c) + 3'i'w' + 34w" 


-ioR((w - c)^\lf - c) 


(17) 
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C is an integration constant. If one now again introduces 


n = (y - ^ = (cxR)"^/3, 


w - c = eatj + \|r = + 


there results 


'l^r 


)"an + 3^0*^ = -i(a^Ti2\lfQ _ c) 


a,Ti + 3>if]_ a + i|fQ bt) + 6i|rQ bi) + bilfQb = 
-i^2abilfQT] + a2^2^J 


(17a) 


Of course, these differential equations still contain all solu- 
tions of equation (7a). We intend to study particularly cpj^ ^ qp 2 

shows for a = 0 at the point Y = Jq regular behavior); therefore, 

we select the one solution of equation (l7a) which behaves at some 

distance from y-,, thus for large (w - c)aE, like i , since 

(w - c)2 

we know from equation (l4) that cp|^ at some distance from Yq is 


given by 


( 


.-c)| 


<iy 


(w - c)' 


Thus we obtain according to equation (l7a) 


= 




and 


Ti + 3\ = -if-p 


0 


(18) 
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is again fully determined by the fact that it is to behave "at 


infinity" like 


2bC 

a3tl‘ 


We now ask for the transformation substitutions for (and cp^); 

meaning thereby the following: In the asymptotically valid repre- 
sentations for 9^ (equation (l4)), we always find the integral 

r dy 
J (w - c)^ 

point y = yQ (w - c = 0). Actually, 


which loses its sense if it is to be extended beyond the 

1 


- is, near the critical 
(w - c)2 

point, replaced by the function i|f. Thus the behavior of ^ (partic- 
ularly ti) in the critical neighborhood is the solely decisive factor. 

If this behavior and therein the magnitude of the integral J' ^ dy 

(extended beyond the critical point) be known, this knowledge is 
equivalent to knowing the transformation substitutions for 9l).* 

The solution characterized by equation (l8) and the boundary 

condition at Infinity reads: 




jtbC 

Sa^Ti 


|®2/3 


( 1 ) 


|(-iaoTl)3/2 




[ +00 


^/3 


(2) 




(19) 


Therein Hankel’ s cylinder functions of the index 2/3 and the 
argument ^(-iaQTj^3/2 appear ^clq = a^/^). The sign of (-iaoTl)^/^ 

is to be taken so that becomes positive for 




(«l)/2 


A closer investigation of equation (19) shows that behaves in the 

entire upper semiplane, and partially, even in the lower one, namely, 

for T] = rei^ within the limits - — < I < at infinity like 
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, if or a is positive. If a is negative^ the upper and 

lover semiplane are interchanged 


lim valid, if 


r — >00 


a^T) 


Hence we infer the important result: 


1 


5i« 13«i 

a < 0 


or 


Titi 


1 


iTt 

6 


a > 0 


( 20 ) 


<iTl = 




2bC 


in 


2bC . 
in 


a > 0 


a < 0 


( 21 ) 


Thus the transformation substitutions for accurate up to the 

magnitudes of the order (aJR)”^^^, are now found for finite values 
of a also; we now know - and that is sufficient - what, according 


to equation (2l), the integral 
to w - c > 0, signifies. 


(w - c)' 


-f extended from w - c < 0 


The formulas (l6) also may he derived once more from equation (lT)j 
to the asymptotic solutions (ll) of equation (7a) correspond the 
integrals 


in (i).(2) 

Ti 2/3 



(19a) 


of the homogeneous equation (l8) (C = O). The problem of finding the 
transformation substitutions of the "asymptotic" integrals (ll) and (l^l-) 
is therewith completed with the required accuracy (except for quantities 

( OR) -1/3). 


of the order 
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Fulfillment of the Bounda,ry Conditions and the 
Stability of the Oscillations Corresponding 
to the Solution System I 

Our considerations so far ha.ve been quite independent of the type 
of profile except for a few limitations concerning the singular points 
which had to be imposed on the profile. In order not to lose our- 
selves in an excessive number of different possibilities, we shall 
further specialize the character of the basic flow. The considera- 
tions, however, have much more general validity. We thus assume that 
the bounding walls are represented by the equations y = +1 and 
y = -1, that, furthermore, the wall y = 1 possesses, with respect to 
the other, a relative velocity in the positive X direction (of the 
magnitude w(+l) - w(-l)) and that the laminar flow adheres to the 
walls (which corresponds to Couette*s test arrangement); finally, we 
assime that in the range -1 < y < +1, that is, in the fluid, Re(w - c 
once and only once is zero. Moreover, we shall presuppose in the 
entire region continuity for w and the derivatives of w and, beyond 
this, make the additional assumption that the functions w, w' , w", 

etc., always are of normal magnitude, that is, that they do not, for 

instance, at certain points, assume a magnitude of the order (olR)^^^. 

Furthermore, for the following calculations, we at first regard 
a as so small and aR as so large that we may put with sufficient 
accuracy 


CP3 


= w - c 


Cpi^ = (w - c) 


dy 


-1 (w - c)‘ 


(l^a) 


The fixing of the lower limit of the integral in obviously does 

not signify a limitation of the generality of our solutions. Rather 
we determine thereby cpj^ as that linear combination of qp^ and cp^^ 

which disappears at the point y = -1. In case w - c should dis- 
appear there also, obviously is replaced by the function 

CP 3 = w - c which now for y = -1 becomes zero. 
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In order to satisfy the boundary conditions first at the wall 
y = -1; we form two aggregates f^, f^ from cpj^ for 

which really cp = cp' =0 for y = -1 


fl = qP4 + 


go(-l)[w(-l) - c] 


w( -] 


cp3 cpl 


I w'(-l) - c 


^2 = ^4 


go(-l)[w(-l) - 0 + |w'(-l) 


r !3_ 

[w(-l) - 


9x(-l) 


n 

9l(-l) 


> ( 22 ) 


Therein we understand from now on by g^ the root ^ -i<xR(w - c) , not 

as in equation (lO), \j-i(v - c), in order to save writing down the 
factor ^oR. 

In order to satisfy the boundary conditions at the other wall as 
well, one must attempt to determine two constants A and B so that 


Af^(+1) + BfgC+l) = 0 


Af^'(+1) + Bf^'C+l) = 0 


The condition for the possibility of such a determination is 


f^(+l) f2(+l) 

f^'(+l) f^’C+l) 


0 


(23) 
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By this condition c or 3 , respectively, is determined if R and a 
are given. Thus it is now a question of solving equation ( 23 ) for c 
and of determining the sign of the imaginary part of 3. Equation ( 23 ) 
forms the perfect analogue to Sommerfeld' s turbulence equation for the 
linear profile. 

From equation (16) we infer 


fji+1) = qpi^(+i) + 


gQ(-l)[w(-l) - 0 - I w*(-l) 
q>j^(+l) + icp2(+l) j 

9 l(-l) J 


q> 3 (+i) 

w(-l) - c 


f3^’(+l) = ipi^'(+l) + 


go(-l)[w(-l) - {] - I w'(-l) 

cpj^’(+l) + iqp^'(+l) 


<P3* 

j^(-l) - ( 






f2(+l) = qpi^(+l) - 


qp^(+l) 


go(-i)|^(-i) - 3 I 


<P 3 (+ 

+ ^ v*(-i) 


1 ) 


- C 


<Pq(-1 ) 


f^u+l) = 94’ (+1) - 


go(-l)|w(-l) - cj + I w'(-l) 


93' (+1) 

w(-l) - c 


^2’ (+1) 
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We insert these values of f2_, ^ 2 ’ ^ 2 * into equation ( 23 ) 

after having made an estimate of the magnitude of the individual terms 
in order to eliminate unnecessary complications of the calculation by 
writing down unessential terms. For this purpose we note that there 
will be^ in general^ either cp2(+l) « 92 ^”^^ 92^”^^' 

This is caused by the factor in the exponent of 9^, cpg in 

equation (ll) if there does not exist the equality 


Re 



c)dy = 


Re / \l-i(v - c) dy 

Uyo 


which we exclude. 

Which one of the two cases will occur cannot be decided before- 
hand; generally^ both are possible and yield both solutions. In the 
case of an obliquely symmetric profile ^ one case gives the solutions 
symmetrical to that of the other. At any rate^ the two possibilities 
behave principally quite analogously and it is therefore sufficient 
to investigate one of the two. Thus we assume 


cp2(+l) « 92(-1) 


that is (compare page 9 )^ the point w = c is to lie nearer to 
w = w(+l) than to w = w(-l). 

Hence it follows that 


1 

is extraordinarily small, thus is very large. Thus there 

9l(-l) 

remain in f]|_ and fq' only the terms which have q 3 q(-l) in the 
denominator; in f 2 and f2' the terms containing 92 elimina,ted. 
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From eqviation (23) we thus obtain 



Even in this form the equation for c is still rather complicated. 

We therefore f\irther simplify eq\aation (25) by cancelling now not only 

quantities of the order of magnitude e but also qxiantities of the 

order (oLR)"^/^. 

For this purpose we determine that gQ(+l) is of the order 

(clR)^/^, thus at first excluding the possibility of w(+l) - c being 
very siTiall, and that furthermore 


cp^'(+l) + icp2'(+l) 


5 w* (+ 1 ) 

^ w(+l) - c 




+ 


So 




Thus we retain only those terms -of equation (25) which are multiplied 
by the factor gQ(+l). 
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Thus the simple result is found 



icp2(+li] <Ph(+l) 


= 0 


or 





-i<xR(w-c)dy 


Pi 


-1 


dy 


(w - c)^ 


= 0 


(26) 


This equation possesses two completely different solution systems 

■>y=+l 


r* I — 

ij 


2 

e -'yo 


c)dy 


= 1 


(I) 


dy 


Li 


= 0 


(II) 


System I represents the perfect ana,logue to the solutions Hopf 
obtained for the linear profile and has discussed in detail elsewhere, 
section Actually it is shown that the oscillations corresponding 
to system I always are of stable character. From 

\| -i(xR(w-c)dy 
e = i 

follows 




2 


( 27 ) 
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where n is a positive (compare page 9) integer. It is easily seen 
that this equation can only be satisfied when ac = p possesses a 
positive imaginary part. Thus the oscillations characterized by eqtia- 
tion (27) are actually damped, the amount of the damping being of the 
order of magnitude w(+l) - c, and therefore by no means need be small. 


5 . The Solution System II and the Conditions for 
Instability of a Profile 

The solutions in the system II are identical with the solutions 
of the Rayleigh equation (8) and satisfy the condition 



(23) 


or (compare the remark to eqioation (l^), page 19) quite generally 


for 



y = +1 


and 


y = -1 

The latter form differs from the first in certain exceptional cases 
which will be discussed later; moreover, equation (28) represents, of 
coxarse, only a first approximation (a = 0). For the solutions of 
equation (28) one must now distinguish four different possibilities: 
Either, (l) equation (28) has solutions with complex c, then the 
profile is always unstable since the conjugate complex value of c 
also always represents a solution; (2) there exist solutions of equa- 
tion (28) with real c. Then we designate, as Prandtl did elsewhere, 
the profile as "capable of oscillations." This can, according to 
eqmtion (2l), occur only if, at the point (w = c), w" = 0, if, 
therefore, the profile either possesses a point of inflection or is 
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composed of linear pieces; (3) real values of c exist which make 
at least the real part of 



zero; or, finally, (4) if none of these three cases occur, eqxaation (28) 
has no solution. In cases (3) and (4), we call the profile "not 
capable of oscillations." We contend that case 1 always results in 
instability, cases (3) and (4) always in stability, case (2) generally 
in instability of the profile taken as a basis. For cases (l) and (4), 
this has already been proved above. In case (3)> we put c = Cj. + icj^ 

with Cj. signifying that real value of c for which the real part of 



disappears. Then we know from section 3 that for 

2"b 

part of the integral becomes jti, for c,- » 

I O I 


Ci ^ 0, 


the imaginary 
, however. 


2b 

I a^l 


«i. 


Thus, for reasons of continuity (compare section 3)^ a 


point > 0 must exist where the imaginary part of the integral (28) 

also disappears. The four solutions of equation (28) thus character- 
ized yield therefore a quantity c with a positive imaginary part, 
thus stable oscillations. 


Case (2) finally requires somewhat more detailed calculations. 
Before performing them we note that to case (2) pertain two types of 
solution for equation (28) which cannot be represented in the form 



If w(+l) = w(-l) a solution of eq\iation (28) is 9 = w - w(+l); in 
fact, here 9 = 0 for y = +1 and y = -1. Furthermore, it 
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may happen that, for instance, w* becomes infinite for w = w(+l). 
Then 

cp = [w - w(+i^ r - — g 

J_1 [w - w(+l^ 

also is a solution of equation (3) which satisfies the boundary con- 
ditions. We shall not treat this case here in more d.etail since it 
will be discussed more thoroughly in part II; however, it must be 
noted as essential that the difference between cases (2) and (3) is 
very large and that it is, for instance, by no means sufficient to 
approximate, according to Rayleigh, a cui^ed profile by a polygon. 

For an illustration of this difference 


Re 




dy 


(w - 


c) 


2 


is represented qualitatively as a function of c in figure 1 where 
the solid curve corresponds to the curved profile, the dashed curves 
to the one consisting of linear pieces. One sees that every break 
causes a new root Re(j) = 0 because J at the point c = '^hreak 

the broken profile varies like i . This corresponds to 

^ ~ ^break 

Rayleigh's well-known theorem that there are as many oscillation roots 
as breaks. Nevertheless, the curved profile does not possess an oscil- 
lation root. After this comment, we revert to our contention that the 
profiles capable of oscillation generally become unstable if the fric- 
tion is taken into consideration. 

For a proof of this instability, we return to eqviation (25) and 
to the more exact solutions in system II. Since we know that c is 

real, except for quantities of the order of magnitude (aR)~^^^, we 
may assume 


cpgC+l) »qP]_(+l) 
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If we furthermore neglect the terms of the order (clR)”^/^ in equa- 
tion (25); we obtain after slight transformations 


m 




1 ( 


w - c, 


1 

go(-l) [w(-l) - 3 ^ 


1 

go(+i) [w(+l) - ^ ^ 


We put further c = Cq + 6 where Cq is- the zero of J, S a small 

-I/2 

quantity of the order (oR) ' . We assume for reasons of simplicity 
a to be positive; then we may on the right side replace c by Cq 

and may develop the left side into a Taylor series in 6. Thus there 
results, if we break off the Taylor series with the second term which 
we presuppose as sufficient approxiraabion 


J(c) = J(cq) + 6 


and from equation (25) because of 


dc 


(=‘"o) 


J(cq) =0; = ^-ioB(w - c) 


(concerning the sign, compare page 9) 


b — \|2aR = 

Ant \ ’ 


1 - i 


1 -I- i 


^'(c=co)'' [co - .(-1)]’/^ [w(-1)-cJ5/2 


Hence follows, because of 


w( -1) > w(+l) - c. 
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(cq is to lie nearer to w(+l)) that the imaginary part of 6 and 
thus also that of c and of |3 has the same sign as 

and that oscillations corresponding to a negative value of — have 

dc 

an unstable character. If therefore our partly linear profile still 

has the property that — < 0 at the point v = c, it is unstable. 

dc 

This condition — < 0, however, is satisfied very frequently, for 
dc 

instance, always when the point w = c lies near one wall (for 
instance, y = +l) and the profile is linear from the point w = c to 
the boundary. 

Summarizing, we conclude; The instability or the stability of a 
profile can be decided for all profiles considered so far by their 
behavior in the case of frictionless fluid. Profiles which are capable 
of undamped oscillations in the latter case and where the friction is 
taken into account become, under certain presuppositions, \mstable. 

The latter profiles must have very special properties as shown above; 
they must, for instance, be partly composed of linear pieces or they 
must have a point of inflection w" = 0. (Compare above.) 

At the same time, however, these profiles of type 2 are the only 
ones still to claim physical interest since they are the only ones 
whose behavior with respect to their stability corresponds approximately 
to Reynolds’ conjectures. Following, we shall show that these profiles, 
in general, really have a critical Reynolds number (with the exception 
of the broken profiles). 


6. The Reynolds Number of the Stability Limit; Numerical 

Calculation on the Parabola Profile 

If, therefore, a profile is prescribed which, for frictionless 
fluid, permits imdamped oscillations and with friction is unstable, 
the question arises, for what minimum value of the Reynolds number does 
instability occur? The simplified equations (25), (26), etc., do not 
suffice for answering this question. We must revert to equation (23) 
and to the forms (ll) and (l4) for the integrals 9^, 

cpj^; however, it is, of coxirse, quite impossible generally, for an 
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arbitrary profile to represent the critical Reynolds number as a 
function of w and of integrals over w; it will only be our task to 
indicate the way by which one arrives at the critical velocity and 
then to perform the calculation on a special example. 

Since in our last calculations a and R had appeared only in 
the combination aR (because we had assumed as small), these 

calculations can yield at best a critical value for oR only, not 
for R alone. Thus we must first investigate the behavior of the 
roots of equation (23) for increasing a^. Of the roots of equa- 
tion (23), only those in the solution system II which satisfy the 
equation cp|^^(+l) = 0 are of- interest. 

Instead of equation (23) we must therefore discuss the equation 



cp = 0 for y = -1, y = +1 


If the profile consists, as in Rayleigh's example, of linear pieces, 
there exists (compare page 2?) always a root of equation (28a) for 
every break and these roots remain in existence for every value of 
o,^. Thus, the broken profile yields no maximum value of and 

therefore cannot ever lead to a critical Reynolds number. 

This is different if (cf. pp. 26 and 27) a solution of equations (28) 
or (28a ) , respectively, with real c is possible for the reason that 
either somewhere in the profile w" =0 or that w(+l) = w(-l), 
cp = w - w(+l) represents a solution of equation (28). These latter 
types of solution always yield a solution of equation (28a) only for 
a very definite value of a^. For w" =0, c is determined by the 
very fact that for w" =0, w is to be w = c; thus the equation (28a) 
defines a quite definite value of however,, for the case 

w(+l) = w(-l) a solution of equation (28a) obviously exists only 
for = 0. 


In 

°It is still presupposed that R and oR are large and a « R. 
Thus critical Reynolds n\imbers will possibly appear if these presup- 
positions are no longer valid; however, the respective Reynolds num- 
bers R would then probably assume values so small that they certainly 
would be of no physical significance. 
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For this type of solution of equations (28) or (28a), which are 
characterized in the limit R = 00 by a very definite value of we 

shall expect that, with the friction taken into consideration, a also 
may vary from its definite value only by small amo\mts. For these 
profiles the appearance of a maximimi value (and in the case w" = 0 
also of a minimum value) for a is very understandable. Thus all 
oscillations, the wave length of which is smaller than a certain 
critical wave length, are in such cases damped for all values of clR. 

After having found an upper limit for. a^, one will attempt to 
determine the approximate magnitude for the lower limit of otR. A 
simple investigation of equation ( 25 ) shows that essential variations 
in the imaginary part of c occur only after the exponent of e in 
the approximate representation (ll) in (pj_(+l), <p 2 (+l) has decreased 

to values of the order of magnitude 1 ; however, if this is the case, 
we very soon reach the critical value (for which the imaginary part 
of c is changed from negative to positive values) as will be shown 
in the numerical example. If we assume that w is essentially linear 
between w = Cq (cq = real part of c) and w(+l) the condition for 
the approximate magnitude of oR reads 

(oR)^/^(w(+l) - c^^/^ 

!= ^ 1 (29) 

w*(+l) 


or 



V* (+1)^/3 

w(+l) . Cp. 


Since in the cases of interest to us w(+l) - c^ will probably be 

small^ we may by assumption form a conclusion as to high critical 
Reynolds numbers. At the same time we note that for a certain value 
of R there will always exist not only a maximum value but also a 
minimum value for a of the unstable oscillations. This follows from 
the fact that we did find a minimum value of oR (not R). 

As numerical example for our general calculations made so far, we 
select the parabola profile because it is physically the most 
interesting one. It is to be classified as "profile capable of 
oscillation" of the type w(+l) = w(-l). 
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Here too we shall consider only the two-dimensional motion, that 
is, not Poiseuille's flow in tubes but the flow prevailing between two 
parallel walls at rest (y = +1, y = -l) under the influence of a 
constant pressure gradient. Thus we put 


w = 1 - y 2 ( 30 ) 


The symmetry of w and v - c permits the deduction that cp must 
he an even function of y.^T Thus we single out, from among the solu- 
tions of equation ( 7 a), two symmetrical particular integrals and attempt 
to satisfy the boundary conditions at one of the walls, for Instance, 
y = -1. Those at the other wall then are fulfilled automatically. 
Obviously we may take simply cp^ those symmetrical integrals. 

For the other we choose 


q5]_(y) cp2(y) 

cp^(o) cp2(o) 


It follows from equation (29a) that for our profile near critical 

velocity c will be sraa,ll of the order (olR)”^"^^; in the following 
calculations we shall thus cancel terms of higher than first order 
in c. Furthermore, we state that will be » qPj^(O) so that 

in the neighborhood of w = 0 and w == c the second symmetrical 
function cp simply is reduced to cp^(y). From equation (l 6 ) then 

follows that we have the two integrals cp^ and cp^ - iqp^ at disposal 

for fulfillment of the boundary conditions for y = -1. Equation (23) 
is therewith transformed into 


17 

If one divides cp into a part even in y and a part odd in y, 
each part of cp by itself must satisfy the differential equation ( 7 a) 
and the boundary conditions because of the symmetry of w - c and w. 
For the general stability investigation of the profile 1 - y^ it is 
therefore sufficient to treat the two cases "cp even" and "odd" 
separately and only these two cases; however, it may easily be seen 
that the assimaption of symmetrical oscillations, that is, "cp odd" does 
not lead to a solution of equation (23) and thus not to unstable 
oscillations. The assumption "cp even" therefore suffices for the 
stability investigation. This is noteworthy insofar as, accordingly, 
all symmetrical oscillations are stable and only unsymmetrical dis- 
turbances unstable. 
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# 


CP3(-1) 

qP3’(-l) 


cp^(-l) - i(p2(-l) 

cpi'(-l) - i(:f^'(-l) 


= 0 


(31) 


In cp^ one must always take y = 0 as lower limit for the occurring 
integrals in order to guarantee the symmetry of qp^. Furthermore, we 
shall develop in 93 only up to magnitudes of the order and in 

the development in (ctR)“^ break-off with the terms of the order 
(cxR)"^. We now write equation (31) in the form 


<Pl'(-l) - icp2'(-l) 93' (- 1) 

9i(-l) - i92(-l) Cp3(-1) 

If one inserts equations (ll) and there results 


r V-icuR(w. 

dvo 


c)dy 


+ 1 


2 I ^ -iciR(w-c)dy 
J - i 




iouRc = - — + 
2c 


(32) 





dy(w - 


c)^ + 



1 




+ a 



( 33 ) 


Since c becomes very small, we assume in first approximation w 
from 0 to c as linear; w ~ 2(y + l). Then we obtain 



NACA TM 1291 


3 ^^ 



(31^) 


If we put 


z . i c3/2(2cdi)l/2 

there arises from equation (33) 

r 

^-( 1 + 1 ) z ^ ^ 3z(l + 1) ^ _ 9 ^ !J2 

^-(l+i)z _ . 2 2c 


dy(w - 


1 + 



+ . . . 



(35) 


This equation is perfectly analogous to equation (26). We are 
interested, above all, in the limiting value R or z, respectively, 
for which the imstable oscillations are transformed into stable ones; 
thus the imaginary part of c is exactly zero. This limiting value z 
will, of coiirse, also be a f\mction of a. Thus we now assume c as 
real and thereby obtain the limiting value of z or R, respectively, 
as a function of a. The minimum value of R on this R(a) curve 
will be denoted as the characteristic Reynolds number for the parabola 
profile. Detailed calculation shows that one obtains by means of the 
form ( 35 ) of the stability equation only the upper part of the curve 
R = R(a) (solid line in fig. 2) which was to be expected according to 
the deliberations of section 6; however, one can calculate the lower 
part of the curve R = R(a) only by using for cp^, approximations 

other (compare equation (l9a)) than the asymptotic formulas (ll). The 
critical Reynolds number denotes just the range where the asymptotic 
formulas cease to be valid. Since this circumstance would lead to very 
complicated nxmerical calculations and since we cannot attach, in 
general, (compare section 7) any essential physical significance to the 
type of instability here characterized, we used rough estimates for 
calculation of the lower part (dashed line in fig. l) of the curve 
R = R(a) which, of course, cannot yield qmntitative results; however, 
the qialitative behavior of the curve is surely reproduced correctly. 
Thus, we conclude from figure 2: 
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1. There exist both a maximuin value of a and a minimum value 
of R for instability. 

2. For a definite value of R there exists a maximum as well as 
a minimiim value of a; Instability prevails within these values, sta- 
bility outside. 


3 . The maximum value of a lies approximately at a = 0.7 

the magnitude of the minimum value of R is of the order 

of 103. A calculation of this minimum value with some degree of 
exactness is not possible from the figure. 


T. Physical Discussion of the Results of Part I 

Let us siommarize once more in detail all results found concerning 
the stability problem. Above all, it became clear in the coirrse of 
the calculation that the problem of the stability of a profile for a 
viscous fluid can generally be decided by treating like Lord Rayleigh 
the limiting case of frictionless fluid (equation (8)). Profiles which 
are unstable then (that is, for R =») remain so for sufficiently 
large finite values of R (section 4) as was-, to be expected beforehand. 
Likewise profiles which, in the frictionless case, are not capable of 
oscillations are found to be stable (section 4) and profiles, which 
according to the investigation by equation (8) permit undamped oscil- 
lations, generally to be \mstable (section 5). This latter case 
obviously is the only one which, physically, signifies something new 
compared to frictionless hydrodynamics; however, it should be emphasized 
that this case, contrary to what one might conclude at first from 
Rayleigh's reports, represents an exceptional case. If one disregards 
the possibilities w(+l) = w(-l), ^'boundary ^ " (section 5), it is 

a necessary condition for the occurrence of this exceptional case that 
somewhere in the fluid w" = 0. The broken profiles consisting of 
linear pieces introduced by Rayleigh belong to those exceptional pro- 
files; however, the only permissible conclusion is that curved profiles 
for the purpose of stability investigation may not be approximated by 
polygons according to Rayleigh (page 27). It is true that one may 
find also for profiles curved everywhere (w" ^ O) (section 5) when 
using the differential equations with friction terms oscillations 
which are damped for every value of aR for which, however, in the 
limit R = CO the amount of the damping like (aR)"l/3 tends toward 
zero; thus one has here also undamped oscillations for R = “, However, 
these oscillation roots are lost if one takes the simplified differ- 
ential equation (8) without friction terms as a basis. Insofar, there- 
fore, this also is not a case of exception to the rule according to 
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which consideration of friction, in cases where the frictionless equa- 
tion (8) permits undamped oscillations, results in an excitation; 
however^ as said above, the possibility of undamped oscillations for 
equation (8) must be regarded as an exceptional case. The parabola 
profile belongs to these exceptional profiles (section 6). If we 
investigate further for the unstable profiles concerning the range of 
values for R within which instability occurs, it is found that 
generally, too, only profiles of the last class may lead to critical 
Reynolds numbers, that is, only profiles which permit without friction 
undamped oscillations; however, among the latter, Rayleigh^ s broken 
profiles nevertheless did not yield (compare also page 30, footnote I6) 
a critical Reynolds number. For Rayleigh’s profiles a minimum value 
of (xR does exist but no maximum value of a; therefore, no minimum 
value of R for the neutral stability either (section 6). Only those 
profiles of the last class for which in the frictionless case only a 
definite value of a leads to undamped oscillations (for instance, 
the types w" = 0 for w = c, w( + l) = w(-l)) result in a maximum 
value of a and a minimum value of oR, thus also in a minimum value 
for R. For a definite value of R there exists therefore for those 
profiles a 'maximum as well as a minimiom value of a for the unstable 
oscillations. All these results are in agreement with the stability 
investigations of hydrodynamic profiles made so far.^® 

The question is now how these mathematical results will manifest 
themselves experimentally. It seems surprising that the stable pro- 
files (for instance, Couette’s^ 9 ) and the unstable ones (for instance, 
Poiseuille’ s) empirically show exactly the same behavior. Above a 
certain Reynolds number turbulence occurs in case of sufficient dis- 
turbances; if the disturbances are made as small as possible, the 
laminar profile may be obtained for arbitrarily high Reynolds numbers. 
Especially the last fact, which has been tested by Ekman (footnote 10 , 
p. 2), on the parabola profile seems to contradict the theory for 
unstable profiles; however, it can easily be seen that this contradic- 
tion is only illusory: The smaller the external disturbances, the 

1 8 

-‘-^Compare the reports quoted in the introduction. 

^^However, compare the interesting investigations of Couette's 
motion concerning its stability against three-dimensional disturbances. 
G. J. Taylor, Stability of a Viscous Liquid Contained between Two 
Rotating Cylinders. Phil. Transact, of the Royal Society London 223 ., 
pages 289-3^3, 1922. 

PO 

This possibility of interpreting Ekma.n’ s tests as a sort of 
starting effect has been pointed out to me by Professor Prandtl. I 
should like here to express my deepest gratitude to him for this and 
many other valuable suggestions. 
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longer it takes (particularly for high Reynolds numbers since the 

excitation there is of the order compare section 5) until 

they noticeably influence the motion. Thus it will always be possible 
to postpone, for flow in tubes, this moment so long that the quantity 
of fluid, the stability of which is dealt with, ha.s already left the 
tube when its instability becomes apparent. The Reynolds number we 

calculated could therefore be tested only on a closed system of tubes 

21 

where the same quantity of fluid always flows. On the other hand, 
the tests by Schiller (footnote 9 , p. 2) which show that below a 
certain Reynolds number only laminar motion exists cannot be included 
at all in stability investigations. The original motion here is not 
laminar; one rather deals with existence or nonexistence of a turbulent 
form of motion. At any rate, one may conclude from all these con- 
siderations only that a solution of the turbulence problem by stability 
considerations alone is absolutely not possible. 

Still, the previous investigations may yield important qualitative 
results for our real purpose, the calculation of the turbulent motion. 

If we interpret the turbulent motion as a cer,tain basic flow with 
superimposed undamped oscillations, we may conclude from our calcula- 
tions that the minimum value R for which this type of motion is 
possible probably also lies at values of the order of magnitude 10^; that 
the wave length of the undamped oscillations lies at 2nh/2, namely a 
at values of the order 1, and that a for a prescribed R is confined 
to certain occasionally very narrow limits; that furthermore these 
oscillations have the character of a wall disturbance as may be con- 
cluded from the smallness of w(+l) - c. These qualitative results 
are quite independent of the special form of the basic flow; however, 
beyond such qualitative criteria the calculations made so far do not 
contribute anything toward the actual solution of the turbulence 
probD.em. 


PART II: THE TURBULENT MOTION 

1. Statement of the Mathematical Problem 


The Reynolds number usually denoted as critical (which is, for 
instance, measured in Schiller's tests and indicates the appearance 

21 

However, the disturbances in stability observed by Ruckes 
(footnote 10, p. 2) for rather small Reynolds numbers are perhaps 
caused by instability according to section 7. This would be quite 
conceivable when the critical Reynolds number according to section 7 
lies below the one for which turbulence (compare part II) is possible 
for the first time. 
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of turbulence in case of sufficiently large disturbances) has no con- 
nection with stability problems and with the laminar flow; it is 
absolutely a characteristic constant of the turbulent motion. Like- 
wise, the Blasius drag law and the well-known conclusion derived from 
it (that the turbulent velocity in the proximity of a wall increases 
with the 1/7 power of the distance from the wall) show clearly that 
the so-called turbulent motion has its own well-defined regularities 
and that it represents a second possible form of motion of the viscous 
fluids. Thus the only way to a solution of the turbulence problem is 
to attempt to eliminate the indefiniteness of the turbulent motion 
and to idealize it until it permits mathematical analysis by Stokes* s 
equations. 

The turbulence problem of hydrodynamics is a problem of energetic, 
not of dynamic stability. There exist two different forms of motion 
of the viscous fluid, each of which has a certain range of values of 
Reynolds numbers within which it is possible. Laminar flow is possible 
from R = 0 to R = » but becomes, however, under certain conditions, 
above a certain value of R dynamically neutrally stable. The tiirbulent 
motion on the other hand exists only above a certain critical value 

pp 

of R and is always energetically more stable than the laminar 
motion. Thus one may in the range of R in which both fonas of motion 
are possible always let the fluid drop from the laminar to the tiur- 
bulent state by means of sufficiently large disturbances. 

In order to make approximate mathematical treatment of the tur- 
bulent motion possible, we again consider the flow between two parallel 
walls and make, first, the following ass\amptions: ^8 The flow is to 
be (a) symmetrical about the X axis with the bounding walls at rest, 

(b) periodic in the X direction with the period^^ 2n/a, and 

(c) periodic with time with the period 2it/|3, and (d) all disturbances 
are to propagate with a speed p/a relative to the X axis, that is 

if the motion is developed into a Fourier series, only products of 
i( pt - ax) are to appear in the exponents of^5 e. 


22compare F. Noether, elsewhere. 

^8rhis statement, too, which represents a simple generalization 
of Sommerf eld' s stability theorem was indicated, for investigation of 
the turbulent motion itself, by F. Noether without further conclusions 
in his paper entitled "Zur Theorie der Turbulenz" (Concerning the 
theory of turbulence), Jahresberichte des deutschen Ph,th. Vereins 23, 
page 138, 191k. 

2^he assumption of a definite a is justified by the result of 
part I that a is confined between certain limits which are, partic- 
ularly in the proximity of the minimum value of R, very narrow. 

need iiuL emjjuaoi^c olic that the actual motions are doubt- 

lessly much more complicated; nevertheless, one may well expect these 
statements to permit qualitative statements concerning the tiirbulence. 
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The Fotirier development of the stream f met ion should therefore 


read 




/ \ /• ^ i(3t-cxx) T -i(Bt-ax) / x 2i(Bt-ax) 

= cpQ(y) + Cp^(y)e + cp^(y)e + Cp2(y)e ' + . . . 


(36) 

The mathematical problem then consists in the determination of the 
odd (according to a) functions cpQ, cp^, and cp 2 (^p, ^ . con- 

jugate to epp, Cp 2 ). At first, the degree of convergence of the 

series ( 36 ) is completely unknown; the question of convergence can be 
decided only after calculation of cpQ, epp . . . If we want to carry 

accuracy so far as to the n"^^ approximation, that is, if we want- to 
calculate cp^ . . . q)j^, we obviously obtain (n + l) simultaneous 

differential equations for the (n + l) unknown functions cp^ . . . 9 ^^. 

Following, we shall need partly the first, partly the second approxi- 
mation. Thus we enter eqmtion (4) with the statement equation (36), 
compare the coefficients of the periodic functions on both sides, and 

break off with the term and thus with 92 ? 92* For 9 q' 

we write w (as in equation ( 6 ) for $'). Thus three simultaneous 
differential eq-uations are produced (the simultaneously obtained con- 
jugate equations need not be written down) 


^J(9p'9p - 9p'<Pp) + 2(92*92 ‘ ^2 ^ 

(9p” - a^9p)(w - - 'Ppw" - 92* (^p" ~ o-^9p) - 


29p’ ( 92 " - ^^ 92 ) + ^p(^2"* ■ " 

“V) 

(’>2" - ‘“^'!>2)(» - !) - »>2 - - 'I’l'I’l"') = 


> ( 37 ) 
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The first of these equations may he integrated twice and yields 

(P^’T 3 _ - Ti'cpi + '^{^2^2 ~ ^2*^2) " ( 37 a) 


C and C]_ signify arbitrary integration constants. 

Since the left side of equation ( 37 a) and w are, according to 
requirement (a), odd in y, C must disappear in our case. 

If we go back from the second to the first approximation, our 
system of equations is reduced to two simultaneous differential equa- 
tions for w and 


Ti'cPi - . ^(w' - C^y) 



(38) 


By way of an interpolation we shall now reflect what replaces eq\ja- 
tion (33) if we do not consider a flow symmetrical about the X axis 
(requirement a), (that is, the flow of a fluid under a pressure 
gradient between two walls at rest), but instead a flow antisymmetrical 
about the X axis (that is, a flow between two walls moved relative to 
each other without pressure gradient as in the Couette case). Require- 
ments (b) and (c) are to be maintained. The statement (36) will then 
no longer be satisfactory since cp-j_, etc., for arbitrary p/a 

are no longer even functions; in order to obtain the entire flow 
pattern in terms of odd functions, we must also include the symmetrical 

oscillations of the form in the formulation for \|f, that 

is, must start with the terms 


% + cpi(y)e 


i( pt-ax) 




Til -ye 


-i( pt-ax) 


Ti(y)e 


-i( pt-ax) 
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As a consequence finally all the terms of the form 
appear in i (elimination of requirement (d)). 

In place of eqioation ( 33 ) there results 


'Pi(y)9i(y)' - cpi(y)'cp;]_(y) + cp^(-y)cp^(-y)' - 


cp3_(-y)'cpi(-y) = - c) 


I ( 39 ) 


(^i" - “^’’i)(" - !) - ‘'"'’>1 = ^ k ’ l " * '‘Si) 


The two equations of the system (38) and ( 39 ); respectively, are of 
simple illustrative significance. 

The second equation is none other than our former stability equa- 
tion (7) which determines the amplitude of the oscillation superimposed 
on a basic flow w and which formed the basis for our investigations 
in part I. The first equation, however, represents the theorem of 
momentum. The left side of this equation essentially indicates the 
momentum transferred on the average by the turbulent vorticity^^, the 
term with w' on the right represents the laminar momentum transfer, 
and the constant C or Cj_y, respectively, is the constant of the 

momentvmi theorem. 

Due to the boundary conditions at the walls = 9^' = 0 . There- 
fore there w' = C or respectively; thus at the walls the lami- 
nar momentum transport surpasses the turbulent one, will gen- 

erally be very large. (Compare the next section.) At the channel 
center, however, that is, in the entire tunnel outside of the immediate 


"^'^e are referring here to the mean momentum in the X direction 
which, for our problem, is transferred in the Y direction. The 
momentum in the X direction equals, on the whole u, the velocity of 
the particle transporting the momentum in the Y direction is v; thus 
the momentum transferred during the unit time uv, on the average uv 
which for the case (36) results in 

uv = -ia(9icp3_' - Ti'Ti) 
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wall proximity, w' is of the order of magnitude 1, thus very small 
compared to C or C]_y, respectively. The turbulent momentum trans- 
port will therefore here completely overbalance the laminar one. 

It corresponds to the structure of the systems (38) and (39) 
that we are able to give immediately a trivial solution of them, 
namely cp^ = 0, w’ = C or, respectively, w' = C]_y, that is, we thus 

revert to the laminar motion. 

Our problem now is to obtain definite results concerning the non- 
trivial solutions of equations (38) and ( 39 )* 


2. The Turbulent Motion in Wall Proximity and 
the law of Resistance 

The most important result concerning the behavior of w in the 
immediate proximity of the walls is the law derived by V. Karman 
(elsewhere) from the Blasius drag law by means of considerations of 

similitude that w in the proximity of a wall increases with ql/T 
(t] representing the distance from the wall). We repeat briefly 
Von Karman^ s train of thought since we are thereby enabled to a general 
visualization of what to expect, even without knowing the Blasius law, 
concerning the behavior of w of the wall. 

As can easily be seen from considerations of similitude, it must 
be possible to represent the shearing stress t acting at a wall 
(that is, the drag) in the form 


T = Kyi ^ f(R) 
h 


where k signifies a certain dimensionless constant. 
If we specially assume a power law there is 



KP 



I 
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From eqviation (to) follows inversely 



The velocity distribution in wall proximity must then be represented 
by an equation of the form 


w 



(Here again (compare equations (l) and ( 6 )) w has been selected 
dimensionless and therefore contains U in the denominator; tj denotes 
the distance from the wall. ) 

We again assume in first approximation a power law (let cj be a 
dimensionless constant) 


w 

W 



€ 


If one now requires the velocity distribution in immediate wall proxi- 
mity to be a function only of t, |a, p, but not of h which is 
physically very plausible, there follows 


1 - I 
1+1 


€ 


(to) 


For 



as corresponds to Blasius* law there results 



In order to understand the physical significance of this result 
we note: w signifies that ^ infinite at the 

boundary, that therefore w infinitely clings to the wall; however, 
it is clear that actually w* at the wall cannot be infinite since 
w’ , on the contrary, essentially denotes the shearing stress at the 
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wall and is therefore, according to eqioa.tion to be equated 

to except for a numerical factor independent of R. 


w 


edge 


R- 


(42) 


w* at the wall is therefore very large for the large values of R 
which are of interest to us. Corresponding to its derivation the 

velocity distribution w ~ will be strictly valid only in the 

limiting case of infinitely large distance from the wall or of friction- 
less fluid (R = ~). These facts can be still more easily comprehended 

if the law w tiVt 

is written in the form tj ~ w'^. From the fact 
that the shearing stress is finite we know that the first term of the 
power development ti(w) must be of the form 7]^w. This term, however, 

is very small, essentially eqixalling the reciprocal value of w' and 
thus being of the order R“^ (compare equation (^2)). 

The meaning of the derived law w tj 1/7 

is thus obviously that 

the series development of t](w) is to start with the terms 


T 

T) = 7^w + 7^w + . . . (lt-3) 

where 7^^ is extraordinarily small and that therefore the first 
term 7j^w for somewhat large values of w may he cancelled compared 
with the second 7yw7. 

According to the explanations above we expect independent of the 
validity of the l/7-law for the basic flow of the turbulent motion 
small curvatures at the centers, and in the wall proximity, clinging 
of the basic flow to the walls. 

For such a profile the investigations of part I do not directly 
apply since there w' , w", etc., had been presupposed as finite; 

however, these investigations can easily be generalized to include 
profiles like the one considered here. (Compare section 5, Page 26.) 
Particularly, the solution of the reduced equation (3) (thus lim R = <») 
with satisfaction of the boundary conditions becomes especially simple 
here; the profile characterized just now belongs, according to sec- 
tion 5, page 26, to those capable of oscillation; a solution of eq\ia- 
tion (28) with real c exists. This is extremely important because 
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it shows that the turbulent profiles are always unstable according to 
section 5^ or, in other words, that it is just the deviation (42) from 
the laminar resistance law which makes an unstable profile and thus 
makes turbulence possible. 

The solution of equation (28) is for a = 0 



By selection of the lower limit of the integral we made cp become 
zero for y = -1; that it becomes zero also for y = +1 due to 
selection c = w(+l) can be seen easily from the following 
transformation 



The integral of the right side becomes at the point w = w(+l) 

infinite of lower order than since w' there ( in the 

w - w(+l) 

limit R = oo) becomes infinite. Thus cp^ = 0 for y = +1. By equa- 
tion (44) we have represented in the limit of frictionless fluid the 
amplitude of the turbulent oscillations and derived from the boundary 
conditions the value for p/a, namely p/a = w(+l). It is, however, 
self-evident that the solution symmetrical to equation (44) 



also satisfies the boundary conditions; thus c = w(-l) is valid 
here. In case of the Couette arrangement we therefore conclude from 
equations (44) and (44a) that two mutually symmetrical oscillation 
systems exist, the wave velocities of which agree respectively with 
the velocities of the two walls (w(+l) and w(-l)). 
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For the symmetrical f lov "between two walls at rest , on the other 
hand, v(+l) = w(-l) = 0. From equations (44) and (44a) we then con- 
clude that every integral of the form 


cp = w 



w^ 


satisfies the "boundary conditions; howWer, from 
that cp is to "be odd there results that we must 

of the integral / — y = 0. Thus 


the requirement (a) 
select as lower limit 


cp = w 




(44b) 


In the case of symmetrical flow there is therefore, particularly 
for p/a 

w(+l) = w(-l) = -^ = 0 (44c) 

a 


In the turbulent basic flow, the type of its singularity at the 
walls is of foremost interest to us; thus for the assimption w ~ q®, 
the exponent e. "We shall attempt to show that from the differential 
equations ( 38 ) and (39) respectively in the limit R = «. at least in 
immediate proximity of the wall such a power law with the exponent 
e = 1/7 actually follows. It is true that the domain of convergence 
of the power series used is not established so that the conclusions, 
as far as they apply to the shape of the profile at some distance from 
the wall, are uncertain. We develop cp^_ and w in the neighborhood 

of T) = 0 in integral and positive powers of q - this is possible for 
any finite value - and then inversely q in integral powers of w. 

Thus we are led directly to the formula (43) for q(w). 

We contend, and this is the most important result we shall need 
later, that cpj_ in first approximation may be represented by a series 

of the form 


cpi = U2q 


2 
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where a 2 , clq . . . are real^ . . . purely imaginary con- 

stants; furthermore, w is of the form 


w 


3tT1 




(^ 5 ) 


This contention may be proved for the differential equations (38) 
directly by expressing qp and w in undetermined coefficients if 
the terms clq, a 2 > and Pq, 3^ are prescribed. Thus we 

will, above all, attempt to determine these terms. First, and 

qp^' for Tj = 0 must be zero because of the boundary conditions; thus 

the series for qp^ starts with (oq = a]_ = O). We can verify 

afterward that, furthermore, the following term is eliminated, 

that is, becomes very small compared tc the other terms. By way of an 
interpolation we shall prove here for this purpose by a single approxi- 
mate integration of the second equation (38) that assumes the 

order of magnitude oJR. For = 0 equation (38) reads 


„l”(v . 1) - i 

Whence follows 



The constant A is here of the same order of magnitude as the left 
side of equation ( 46 ) at the center of the tunnel, thus almost of the 
order of magnitude 1. (Compare part I, section 2.) The term cp^'" 

at the edge is therefore of the order oR due to the boundary condi- 
tions. The same is valid for a^. 


Thus we shall meanwhile assume as small and later attempt to 

justify that assumption. Of the constants 3 q, 3^_, the first, 3 q, 

equals zero because of the requirement (a), section 1. 
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The constants and P-j. are, at first, arbitrary^' and there 

is no possibility of deriving them from the solution of the differential 
equations ( 38 ) and (39) in the proximity of the wall. This possibility 
would .arise only if we should succeed in continuing the solution (45) 
analytically up to the other wall; however, this is an extremely com- 
plicated mathematical problem if only for the reason that, as will be 
seen, the simplified equations ( 38 ) and (39) are not sufficient for 
determining and w at the center of the tunnel. Although we 

must therefore forego the solution of this problem, we may still expect 
to obtain, by merely developing and w in the proximity of one 

wall with undetermined coefficients (X 2 , 3^^, those qualitative char- 

acteristics of w and cp]_ in wall proximity which are, according to 
experience, quite independent of the behavior of the fluid at the 
tiuinel center as, for instance, the law w ~ 

We enter eqmtion ( 38 ) with the statement 


<Pl = + a^vp + . . . 

w = 33 _ti + + . . . 


replacing the second equation by (46). We therefore again assme a 
as very small which here only signifies (compare part I, section 2) 
that the wave length of the oscillations is to be large compared to 
the boundary- layer thickness; moreover, we put, according to eqxia- 

tion (44c) 

i = 0 

a 

For the first eq\iation ( 38 ) we write furthermore 


-ioR^Cp^' cp^ - cp^'cp^) = 2aR^(p3_^'qP3_j. - 9ii<Pir') = w' - C^^y 

Therein cp2.r <i®notes the real, cp^^^ the imaginary part of cp^. 

^"^We shall assume 02 as real. This does not imply a limitation 
of generality since cp^ is determined only up to a factor of the form 

e^^ as the initial point of the time coordinate in equation (36) may 
be chosen arbitrarily. » 



From. eqmt ions (46) and (4 t) now follow the recursion formulas 


n(n - l)(n - 2)o^ = -ioJR 



s(n 


2s) °4i-l-sPs-l 


(^) 


nP^ = 2oR s(n 




m 


in addition 

^1 ^l^edge 
23g = 

Therein ag^ denotes the real, the imaginary part of ag. 

From equation (48) there follows first 

0-4 = 0 

From eqxoation (^9) there then results 


Pj = Plj = 3, = Sg - 0 


The tenn ^2 ^7 also he approximately equated to zero. 

From equation (4-9) there follows 
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For very small the term ^2^^ therefore to be neglected com- 
pared to the first term for larger t], however, the higher terms 

3^^”^ etc., are completely predominant. 

Let us thus assume also 32 = 0 and thus calculate the higher 
terms of the series for 3i and w. 

There follows 


023 i 2 

ote = -ioR ; ; = Oy = 0 ; ao = -(oR) ^ 

^ 3x4x5 ° ° 3 x5x6xTx8 

°9 = “-lO = 


= -(clR) 3 ( 




° 23 i- 


2x7x9x10x11 3x5x6x8x9x10x11/' 

g-iP = ® 


a ^3 

S7 ' ' %o ■ - ^12 = «; 


3i 3 = (oR) 




02^1 


b X 5 X 6 X 8 x 11 X 13 X l4 7 x 10 x 11 x 


rJ 


The representation ( 45 ) for w we contended is therefore proved and 
it may easily be shown too that of the further terms only in every 
sixth term has 3 a finite value. 

Hence follows for the representation of as power series in w 


T] = 7^w + 7 .j,w‘ + 


13 + 

' • • • 


77 - (or) 


2 

703 i 


M50) 


r = (art)** / X 3 X 3 X 17 

^3 \7 X 10 X 10 X 11 X 13 X 3^13 




5x5x6x8xiix13xi4x 3i^^y 
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The terms to y^, 7g to y-.^, 7-,i. etc., all eqioal 

28 

zero. 


The development (50) now actmlly completely agrees with equa- 
tion (^3) and we seem thus to arrive, even without knowledge of the 

constants 02 and to the law tj w"^ semiempirically derived 

by von Kerman. The coefficients 7^^ and yrj, however, cannot be 

calculated. Inversely, we may perhaps conclude from the empirical 
findings for the coefficients and 02 that 7y is of the order 

of magnitude 1, y^^ of the order thus 02 ~ 

Subsequently, we thus also confirm our former assertion a^ij 

Raising the question of what order of magnitude are the values of w 
for which the third term in equation (50) is small compared to the 

second for which therefore the w^ profile actually is valid, one 
finds w thus '^R" 1/8. 

Accordingly, the profile w ~ 

follows from the differential equations (38) only qualitatively at 
first. No information about the fact that the 1/7 profile has been 
observed almost up to the tuimel center is given in our calculations; 
however, this was not to be expected since the other constants entering 
the law also depend on the behavior of the fluid at the opposite wall. 


As an interpolation, we shall once more briefly summarize what 
factors we have neglected in deriving equation (50) from equations (W) 
and (^9) and attempt thereby to determine within what acc\xracy the 
conclusions drawn from (50) are correct. First, we used system (38) 
instead of (3T)> thus cancelled magnitudes of the order <P2/w. Further- 


more, we eqxiated ag = 0, p/a = 0, ^2 ~ ^ therewith neglected 

aoTj o 02^ 3l 

magnitudes of the order — =— , — , , and respectively. 

02 aw Pi 

The accuracy of our calculations will be determined by the largest of 
the terms here neglected. Simple considerations of the order of magni- 
tude, not executed here, make it probable that of these terms 92/'^ 


is the largest but that this term, too, goes toward zero with R — 


This power series ti(w) may, of course, also be derived directly 
from equations (46) and (4 t) without the detoxir over the series of w(ii) 
if w is introduced as independent variable; however, the calculations 
required for this purpose are somewhat more complicated. 
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Selection of a sufficiently large value for R will therefore make 
it possible to carry the accxiracy of the results derived from equa- 
tions (48), (4-9), and ( 50 ) arbitrarily far. 

As to Blasius' law of resistance, it can, of course, be derived 
inversely according to the method described above from the law q ~ w"^ 
by means of consideration of similitude if one assumes, as we did, that 
the behavior of w in the proximity of the wall is independent of the 
tunnel width; however, for the reasons stated above (impossibility of 
analytical continiiation) we must leave the question unanswered whether 
this latter - physically very plausible - assimiption also follows from 
the differential equations ( 38 ) and ( 39 )^ respectively. 

We are, however, able to draw a noteworthy direct conclusion con- 
cerning the law of resistance from eqviations ( 38 ) and (39) hy means of 
consideration of similitude. In the tunnel, except for immediate wall 
proximity and the point y = 0 (compare below equation (66)) one may 
write instead of the first equation ( 38 ) because of the magnitude of 
(compare pages 4-1 and 4-2) 


ioR(cp3^’(P3^ - cp^'cpi) = C^y (51) 

Since the amplitude cannot go toward infinity with R — ^00 - this 

would render all our calculations devoid of physical sense - there 
follows that is at most of the order of magnitude oR, that there- 

fore the exponent | of equation (4o) must be <1 (which in a certain 
manner also can be seen from equation (4-1)). Hence follows (compare 
eq-uations (42) and (4o)) that the law of resistance t = const. 
usually assumed in hydraulics represents an upper limit for all imagi- 
nable laws of resistance of turbulence which is independent of the wall 
characteristics. One may conclude as an assumption that the law 

7/4 

T ~ U ' is valid only for smooth walls - it was for those only that 
we obtained q ~ w^^^ - that the law of resistance for rough walls, 

» pQ 

however, more and more approaches the qioadratic law. For rough walls 
the amplitude will be independent of R and of the magnitude of 

the wall distvirbances; moreover, for rough walls the boxmdary conditions 
will no longer cause cp^ to be real in first approximation as corre- 
sponds to equation (44). 


^Compare the more exact investigations by Von Kerman, elsewhere, 
and the experimental investigations by Schiller, same periodical 3> 
page 2, 1923 . 
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Nothing is changed in the conclusions of this paragraph if the 
equations (39) are taken as a basis instead of the differential equa- 
t ions ( 38 ) • 

3 . The Turbulent Motion Outside of the Immediate 
Proximity of the Wall 

It is essential for the motion at the tunnel center that here 

is composed of those two integrals of (7a) which appear in case of 
frictionless fluid, thus for equation (8). (Compare part I, section 2.) 
The most important characteristic of qp^ following from this fact is 
that it satisfies - except for magnitudes of the order q >2 and 
(otR)"^ - the condition 


qpj^’cp-j^ - qp-j^'cp^ = Const. 


(52) 


This results, according to Abel's theorem, from the fact that, 
except for magnitudes of the order cpg (oR)"^, and qp^^^ 

(the real and imaginary part of qp^^) are solutions of the differential 
equation (8). 

Hence it can be concluded that the equations ( 38 ) and (39) are not 
sufficient for establishment of the motion over the entire tvmnel width 
but that we have to go back to equation (37) and to the system of equa- 
tions which corresponds to it for Couette’s case. 

This, in general, involves a complication of the mathematical 
problem. Only in Couette's case may bhe problem be solved comparatively 
easily because the first eqxiation (39)^ except for magnitudes of the 

order w'/C, thus (oli)"^/^ and qpg^, compare equation (37)^ agrees 

with equation (52). Whereas, therefore, eqioation (52) in consequence 
of its derivation from Abel' s theorem is' correct only up to magnitudes 
of order qp 2 , in Couette's case equation (39) sho^lld still be valid up to 

m ag n itudes of the order qp 2 ^» This requirement is satisfied if we put 
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This equation is therefore to be regarded as solution for Couette* s 
case, of the differential equation we took as a basis. 

According to equation ( 53 ) it would follow for from ( 37 ) 


Tl 9i “ 9]_ “ 0 


(5i^) 


The system applying to Couette' s case is not equation ( 37 ) but a more 
complicated one which we are not going to write down. We do, however, 
state about it that it leads, like (37)> for cp2 = 0 to the solu- 
tion (5^) and thus to the result 


= ae^y + be" 7 y (55) 


Here, a, b, and 7 are any complex constants. For w then follows 
from the second equation ( 39 ) or, respectively, from its reduced 
form ( 8 ) 


w - c = a^^e^^^ + b^e^^^ (56) 


Since at one of the walls there should be w - c = 0, and since, on 
the other hand, w should be odd in the neighborhood of y = 0, it 
follows that w, simply by the vanishing of and a suitable 

Increase of and b^^, must degenerate to a linear profile. 

Thus, we obtain the important result that for Couette' s case the 
basic profile w of the turbulent motion takes an essentially linear 
course over the entire tunnel width - however, strongly deviating from 
the laminar profile, it will be much flatter than the laminar one - 

[that, however, (compare II, section 2) at the edge it clings again 
like to the wall0 . 

We shall now turn to the more complicated case of a flow between 
two walls at rest, thus exactly to the system ( 37 ). For a solution 
we must naturally be content with rough approximations. First, we 
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can cancel in equation ( 37 ) the right sides of all three equations, 
namely the friction terms; this is fully justified by the considera- 
tions of Part I, section 2. Then we equate p/a = 0 (compare 
equation (44c)). 

We thus obtain for CP]^ in the place of the second equation (37) 


cpj_"w - w"cp^ - a^wcp^ - Cp2'(9i" - - 292(91"’ - + 2f^' (92" - 

^^ 92 ) + ^i('P2'" ' ° ^57) 


If we develop 9 ^ as solution of the equation (57) in powers of 
on one hand and powers of 92 on the other, and if we further note 
that 9 i is to be odd (compare (44b)) and write 9 ^ = 9 ^^ + 9 ^^ there 
results with only the linear terms taken into consideration 


9 


10 




(53) 


9ll = aw 


r >r «(•■ 


So * ‘'2''Plo" * S'»lo’ 


2'Pio> 2" - ’’lo’’2'") 


(59) 


Naturally 9 ^ is herein not fully determined - the constant factor a 
assumed as real which does not signify a limitation remains imdetermined. 

If we substitute this value of 9 ^ into the simplified first 
equation ( 37 ), namely 


9l9^" - 9i9i" = Const. 


(60) 
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we obtain, with (P 2 ^ denoting the imaginary part of the result 


'P2i''Plo" 


2in0 


- 2(p^o''P2i" " 


10'^2i 


/Jo v2- 


Const . 


( 61 ) 


Now, however, as follows from the third equation (37) and from the 
fact that cpi is real in first approximation, satisfies the 
equation 


qp2j^"w - cp2iw" = 0 

thus 

”21 I -=’>10 

If we substitute this value of cp 2 i into equation ( 63 ) and if we 

further consider that for y = 0 the left side of equation ( 61 ) and 
therewith the constant on the right side is zero (this signifies for 
the constant of the right side of equation ( 60 ) only that it is in 

o 

first approximation zero, that is, small of the order cp]^<P 2 or 
a^(Pl<P 2 ; respectively, or a^(P 2 f we obtain 


(62) 


(63) 


'^lo'Pio ■ "PlO 'PlO 


= 0 


(6k) 


which fully agrees with (5^). 

This equation, it is true, becomes, like equation (5^)> trivial 
in the neighborhood of the point y = 0; it is there fulfilled 
identically since (p is an odd function of y. Thus it cannot permit 
there a determination of w. This leads for the symmetrical profile 
(64) to a remarkable discontinuity at the point y = 0. (For the odd 
profile such a discontinuity cannot be seen from the differential 
equations.) If one integrates (37a) one obtains, as shown above, 
after a single integration the equation 

2oE|$^"q)j^ - + 2(^2"^2 ~ ^ 2 "^^ + . . = w" - C ( 65 ) 
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where C (compare pp. 1<-1 and h2) , Blasius* law of resistance being valid, 
is of the order of magnitude (aR)^/^, thus at any rate very large . 

The left side of equation (65) disappears, however, with and q>2 

(which, as we know, are odd functions of y) at the point y = 0. Thus 


w 


y=0 


= c 


( 66 ) 


must he valid there. This signifies that Wy_Q" is very large 

(~(oR)3/^) and that therefore w at the point y = 0 shows a sharp 
hreak^*^ (radius of cvirvature ~(olR)" 3/^^. At a small distance from 
this point the coirrse of w must, according to equation (64), again 
he essentially linear. 


We obtain the result: For the flow between two wa^lls at rest as 

well - and siorely this may be applied also to the flow in the tube - 
the profile is linear approximately over the entire tunnel width; at 
the center, however, it shows a sharp break ( it clings to the walls 
with the y^/' law). (Compare figure 3.) 


The physical cause of the sharp break is the fact that the gradient 
of the turbulent momentum transfer for y = 0 disappears for reasons 
of symmetry and that therefore, because the gradient of the entire 
momentum transfer over the tunnel width is constant, the gradient of 
the laminar momentum transfer, that is w", must be very large there. 


4. Final Remarks and Summary of the Physical Results 

Our investigations still show two important gaps. First, they do 
not yield the transition from the t)V 7 profile to the linear profile 
valid in the center part. Second, they are limited to large values 
of R and thus do not yield the minimum value of R, either, if such 
a minimum value exists for which the turbulent motion is still possible. 
The first of these two gaps is most difficult to fill in (compare 
page 48); we cannot even indicate a method which would satisfactorily 


^^Professor Prandtl was so kind as to point out this break to me 
on the basis of empirical material. The break seems less sharp empiri- 
cally than according to calculation results, which is easily explained by 
the fact that the assumption (a), page 38 ^ concerning the symmetry of the 
vortices nnd disturbances also does not exactly correspond to actual 
conditions. 
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solve this particular problem. One may attempt to piece the two 
approximations together that come from the wall and from the tunnel 
center. This would have to be done by means of the condition that 
at the respective junction and w' are to 

be continuous; however, the convergence of the developments (U 5 ) and 
( 50 ) is hardly siofficient thus to guarantee a somewhat defined approxi- 
mation. At any rate the ultimate result, the profile w, is still to 
a great deal dependent on the type of joining the two approximations. 
Finally, it must be regarded as dubious whether such an exact carrying 
out of the formulation (page 39) would yield essentially new physical 
results in agreement with experience since these statements certainly 
represent a very strong idealization of actual conditions. 

In contrast, filling in of the second gap does not offer any 
basic difficulties whatsoever; all necessary expedients are contained 
in Part I and once the profile w is completely known, the methods 
described in Part I are, on principle, sufficient to calculate according 
to Part I, section 6, the minimim value of R for which the turbulent 
motion is possible. One could, for Instance, calculate the critical 
Reynolds number for a profile obtained, according to the method 
mentioned above, by piecing together the two approximations, or one 
could base this investigation on the empirically observed profile and 
thus calculate the Reynolds number in a semiempirlcal manner. In any 
case one will - the investigations in Part I made this probable and 
direct calculations, here not reproduced, confirmed it - arrive at 
the same order of magnitude of the critical Reynolds nixmber, namely 
R ~ 103. The exact value of R will, it is true, still be too 
dependent on the manner in which the profile was obtained to permit 
comparisons with experience. For that reason we did not perform here 
such a calculation of R. 

Let us finally summarize what may be concluded as physical result 
from our investigations concerning the turbulence problem. In Part I 
we recognize that the laminar motion and its stability condition are 
not of essential significance for the turbulence problem and the 
critical Reynolds nimber. In Part II, however, we investigated the 
turbulent motion itself and may hence give a few data on the turbulent 
state of motion. In general, the velocity distribution over the 
entire tunnel is of the simplest type; it is - according to the test 
conditions - linear or constant (section 3)* At the center there is, 
for symmetrical flow between two walls at rest, a sha.rp break; at 
the walls the flow clings, for the ^ 1/7 

profile, to the walls 

(section 2). The calculations do not disclose anything about the fact 
that the 1/7 profile is valid until far into the tunnel interior. The 
turbulent oscillations are for Couette's case almost harmonic in the 
interior of the tunnel (section 3^ equation (53) )j in the proximity 
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of the walls all oscillations will occur. The velocity of the waves 
agrees with the wall velocity (section 2, eqviations (44) - (44c)); for 
Couette's case there exist two groups of turbulent oscillations, one 
of which agrees, with respect to its velocity of propagation, with one 
of the walls, whereas the other group possesses the velocity of the 
other wall. Thus the t\irbulent disturbances show, superficially, the 
character of a wall disturbance. It must, however, be emphasized that 
these distiarbances are capable of existence as free oscillations, inde- 
pendently of wall roughness and similar influences. The amplitude of 
the turbulent waves considerably increases toward the walls (this 
follows from equation (44), section 2) and goes toward zero only 
directly at the wall. 

The wave length of the occurring oscillations (Part I, section 8) 
is, with respect to order of magnitude, equal to (rather somewhat 
larger than) the tunnel width. The minimiun value of the Reynolds num- 
ber (Part I, section 8) for which turbulence is still possible, lies - 
with respect to order of magnitude - near 103. From the profile 

Blasius' T u"^/^ seems to result, under certain presuppositions, as 
the law of resistance for smooth walls. For rough walls it probably 
approaches (section 2) the hydraulic law t ~ u^. The purpose of the 
present report was not so much to establish these regularities, to a 
great part known before, as it was to prove that all results obtained 
so far (seemingly partly contradicting each other) can be uniformly 
described mathematically with the aid of simple basic assumptions. 


I wish to express here my deepest gratitude to my revered teacher. 
Professor Sommerfeld, for suggesting this report and for frequent 
assistance. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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